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The temperature and magnetic-field (H) dependences of thermal conductivity (ft) of 
Bi2Sr2CaCu208+i (Bi2212) are systematically measured for a broad doping range by using both 
pure Bi2212 single crystals with tuned oxygen contents and Bi2Sr2Cai_ a ;Dy a ;Cu208+i (Dy-Bi2212) 
single crystals with different Dy contents x. In the underdoped samples, the quasiparticle (QP) peak 
below T c is strongly suppressed, indicating strong QP scattering by impurities or oxygen defects, 
whereas the phonon conductivity is enhanced in moderately Dy-doped samples and a phonon peak 
at 10 K is observed for the first time in Bi2212 system, which means Dy 3+ ions not only introduce 
the impurities or point defects but also stabilize the crystal lattice. The subkelvin data show that the 
QP heat conductivity gradually decreases upon lowering the hole doping level. The magnetic-field 
dependence of k at temperature above 5 K is mainly due to the QP scattering off vortices. While 
the underdoped pure Bi2212 show very weak field dependence of k, the Dy-doped samples present 
an additional "dip" -like term of k(H) at low field, which is discussed to be related to the phonon 
scattering by free spins of Dy 3+ ions. For non-superconducting Dy-Bi2212 samples with x ~ 0.50, 
an interesting "plateau" feature shows up in the low-T k(H) isotherms with characteristic field at 
1 - 2 T, for which we discuss the possible revlevance of magnon excitations. 

PACS numbers: 74.25.Fy, 74.72.Hs 



I. INTRODUCTION 

Low-temperature heat transport has been extensively 
studied for high-T c cuprates, since it not only pro- 
vides the most straightforward information of quasi- 
particle (QP) transport properties and the peculiar 
electronic state, but also shows rich physics like the 
magnon heat transport of low-dimensional spin system, 
the phonon heat transport and its interaction with pecu- 
liar spin/charge order. 

It was predicted that the heat transport of nodal 
QPs in a d-wave superconductor presents a "universal" 
behavior ^i 3 . that is, the residual conductivity in the 
zero-T limit kq/T is independent of QP scattering rate. 
Physically, it results from the exact compensation be- 
tween the decrease of scattering time and the increase 
of zero-energy density of states with increasing impu- 
rity density. In the standard self-consistent T-matrix ap- 
proximation (SCTMA) theoryji when the impurity band- 
width 7 is in the universal limit ksT <C 7 <C Ao (Ao is 
the maximum gap), kq/T is expressed as 2 * 3 - 

«o _ fc|n fvp_ + V2_\ ^ fc| nvp_ ,^ 
T 3fr c \ V2 vp J 3h c V2 

with n the number of Cu02 planes per unit cell, c the o 
axis lattice constant, and vf (v 2 ) the QP velocity normal 
(tangential) to the Fermi surface at the gap node. 

The early measurements on nearly optimally- 
doped YBa 2 Cu 3 O y (YBCO)^ and Bi 2 Sr 2 CaCu 2 8+5 



(Bi2212)£i£ with various impurity scattering strengths 
were found to be supportive of the "universal" behavior. 
If valid, the universal nature is very useful because 
one can obtain Ao from the bulk measurement of k.— ^ 
Sutherland et alJ analyzed their residual conductivity 
data of YBCO in the framework of the "universal" 
thermal conductivity and obtained the gap size Ao 
from Ko/T; based on this analysis, it was argued that 
the Ao values coincide with those obtained from other 
spectroscopic measurements. Together with the similar 
analysis on the k data of Tl 2 Ba 2 Cu0 6+( 5 (T12201) in 
a wide overdoping ranged a pure and simple <i-wave 
superconducting state throughout the phase diagram 
was proposed. However, as more detailed and accurate 
data are accumulated, it has become obvious that there 
are cases that do not seem to obey the "universal" be- 
havior, including the increase of kq/T with doping level 
confirmed in La 2 ^Sr :E Cu04 (LSCO), 7 ' 9 YBCO^iii 
and Bi 2 Sr 2 _ 2; La 2 ;Cu06+5 (BSLCO)r^ and the varnish- 
ing residual thermal conductivity in an underdoped 
YBa 2 Cu40g (Y124)»i 3 - A possible explanation for these 
discrepancies was suggested based on a systematic study 
of Zn-doping effect on the residual thermal conductivity 
of several typical cuprate families^ By using exception- 
ally high quality single-crystal samples and achieving 
very small experimental errors, it was found that the 
"universal" picture is not a precise description of the 
low-T QP transport properties^ whereas the electronic 
inhomogeneity was found to be the main reason for 
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the failure of the classical d-wave theory^ It should 
be mentioned that there is a controversy over the data 
analysis of the very-low-T thermal conductivity that is, 
whether the commonly observed T-dependence of k/T 
data below 1 K weaker than T 2 is due to the phonon 
specular reflection off the sample surface or simply 
because the phonon boundary scattering can only be 
achieved at temperature lower than ~ 130 mK4^ ' n i 14 
The key question is whether the phonons can be scat- 
tered by something else except for the boundary at very 
low temperatures. An indirect but helpful result for 
this issue was recently obtained in another novel oxide 
material GdBaC^Os+z (GBCO); the phonon scattering 
by free spins in this material was found to be very active 
in subkelvin temperatures, which makes the boundary 
scattering limit not achievable at temperature as low as 
0.3 

The magnetic-field dependence of low-T thermal con- 
ductivity is very informative for the QP transport prop- 
erties. The pioneering work in nearly optimally-doped 
Bi2212 discovered a striking "plateau" behavior of k(H) 
isotherms in the temperature range down to 4 Kiii It 
was originally discussed that the "plateau" behavior is 
due to some kind of field-induced phase transition ; 17 i 18 
which, however, was questioned by a subsequent experi- 
mental observation that the k value in the "plateau" de- 
pends on the history of applying the magnetic fields It 
was then proposed that this behavior can be well under- 
stood by considering the competition between the QP 
scattering by vortices, which is predominant at higher 
temperature) 17 i 20 ' 21 i 22 and the field-induced QP excita- 
tions, which is known as the "Volovik" effect that is 
more important at low temperature i 21 ^ 23 This expla- 
nation was supported by the n(H) measurements at 
subkelvin temperatures ! 24 ! 25 Furthermore, the in-plane 
anisotropic heat transport measurements done on Bi2212 
single crystals discovered that the "plateau" behavior 
is only present for heat current along the b axis, 26 
probably due to the distortion of d x i_ y i pairing sym- 
metry. Subsequently, more efforts were devoted to 
the studying of k{H) behaviors at subkelvin temper- 
ature for cuprate superconductors with different dop- 
ing levelsi 10 i 12 ' 27 i 28 One interesting finding is the field- 
induced suppression of the very-low-T heat conductivity 
in the underdoped La2_ :r Sr. r Cu04 (LSCO), whereas the 
magnetic field naturally enhances the thermal conduc- 
tivity of optimally-doped and overdoped LSCO as ex- 
pected from the "Volovik" effecti 27 i 28 The unusual be- 
havior in the underdoped LSCO can be understood as 
the QP localization associated with the magnetic-field- 
induced spin or charge order j 27 i 29 i 30 ' 31 i 32 i 33 Similar k(H) 
results have been reported for weakly doped YBCO and 
BSLCOi 10 ! 12 Note that in all these previous studies, the 
field dependence of thermal conductivity was attributed 
only to the electron transport and the contribution of 
phonon heat conductivity was assumed to be indepen- 
dent of magnetic field. However, it was recently found 
that in an undoped compound Pri.3Lao.7Cu04 (PLCO), 



dilute (~ 1%) free 1/2-spins (which are presumably cre- 
ated by a small oxygen nonstoichiomctry) can strongly 
scatter phonons and produce a pronounced "dip" behav- 
ior in the k(H) isotherms at low temperatures^ 4 -. Such 
an observation implies that the phonon heat transport 
could be sensitive to magnetic field in high-T c cuprates, 
which are known to commonly possess some local mag- 
netic momentS i 35 i 36 ' 37 i 38 ' 39 i 40 ' 41 Apparently, correctly un- 
derstanding the phonon heat transport property becomes 
crucial, not only for its own research interests but also 
for its role in explaining the magnetic field dependence 
of low-T thermal conductivity and the issue of boundary 
scattering limit, which are fundamental for capturing the 
QP transport properties. 

In this work, we study the low-T thermal conductivity 
of high-quality Bi2212 single crystals at various dopings 
to probe both the QP heat transport and the phonon 
heat transport properties. It is worth mentioning that 
the doping dependence of the heat transport has been 
seldom studied for Bi2212. 14 Besides confirming the de- 
crease of QP heat transport with decreasing doping by 
either tuning the oxygen content or doping Dy, we obtain 
several novel results: (i) Dy substitution for Ca stabi- 
lizes the crystal lattice and leads to the improvement of 
the phonon heat transport; (ii) the low-T phonon peak 
in k(T) curves is observed in appropriately Dy-doped 
samples, which is the first time for Bi2212 system; (iii) 
Dy doping introduces some magnetic-field dependence of 
phonon heat transport and a "plateau" phenomenon in 
k(H) isotherms is found in the non-superconducting Dy- 
doped samples. 



II. EXPERIMENTS 

High-quality pure Bi2Sr2CaCu2 0s+5 and Dy-doped 
[Bi2Sr2Cai_ I Dy I Cu208+a] single crystals are grown by 
the floating zone (FZ) metho d 22 ! 26 . The pure crystals 
are carefully annealed at 400 - 800 °C and in different 
atmospheres to tune the oxygen content, e.g., the opti- 
mal doping with T c = 95 K is obtained by annealing at 
800 °C in air. Most of the Dy-doped samples, if not spe- 
cially mentioned, are annealed at the same condition as 
that of the optimal doping. We label the pure samples 
as UD70, OP95 and OD70, etc., for showing their doping 
regimes and T c values, and Dy-doped samples as Dy81, 
Dy45, etc., by their T c values (the Dy doping effectively 
reduces the hole concentration and brings the samples 
into the underdoped regime). The actual Dy concentra- 
tion x is determined by the inductively-coupled plasma 
atomic-emission spectroscopy (ICP-AES). 

The thermal conductivity k is measured along the a 
axis, perpendicular to the direction of the incommen- 
surate modulation structure. The typical size of the 
sample is 2.5 x 0.8 x 0.03 mm 3 , where the longest and 
shortest dimensions are along the a and c axis, respec- 
tively. The measurements are done by several different 
processes: using a conventional steady-state technique by 
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FIG. 1: (color online) (a) In-plane resistivity of several pure 
and Dy-doped Bi2212 single crystals, (b) Renormalized Hall 
coefficient of these crystals. LSCO data (solid lines) from 
Ref. are included for comparison, (c) T c vs p for several 
pure and Dy-doped Bi2212 single crystals. Doping depen- 
dence of T c does not follow the empirical formula T c /T™ ax = 
1 — 82.6(p — 0.16) 2 , shown by the dashed line. 
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using Chromcl-Constantan thermocouple in a 4 He cryo- 
stat for k(T) data at "high" temperatures (2 - 150 K) 
and for k(H) data above 40 K; using a "one heater, two 
thermometers" technique with Cernox sensors in the 4 He 
cryostat for k(H) data at 5-30 K; using a "one heater, 
two thermometers" technique with RuC>2 sensors in a 
dilution refrigerator or a 3 He refrigerator for k(T) data 
at subkelvin temperatures and for n(H) data below 10 
K. All the measurement techniques for k have been de- 
scribed in our previous publications^ ' 10 ' 14 ' 26 ' 27 ' 34 The er- 
ror and irreproducibility of thermal conductivity data are 
typically less than 10 %. 

Magnetization measurements are carried out using a 
Quantum Design SQUID magnetometer. 



III. RESULTS AND DISCUSSION 
A. Electrical transport and carrier concentration 

Figure 1(a) shows the in- plane resistivity of several 
Dy-doped Bi2212 single crystals. The optimally-doped 
sample (without Dy doping) shows a "negative" resid- 
ual resistivity^ suggesting very high cleanliness of our 
crystals. Upon doping Dy the superconductivity is grad- 
ually suppressed and the increase of residual resistivity 
is rather weak, thus the main reason for the suppression 
of the superconductivity is the reduction of hole concen- 
tration when substituting Dy 3+ ions for Ca 2+ ions. 

It is difficult to determine the carrier concentration in 
Bi2212 because the hole concentration per Cu, p, cannot 



FIG. 2: (color online) Temperature dependence of thermal 
conductivity of pure and Dy-doped Bi2212 single crystals. 
The Dy contents x for Dy81, Dy45 and DyO samples are 0.17, 
0.30 and 0.53, respectively. 



be chemically determined like in LSCO. Ando et al. have 
shown that the magnitude of renormalized Hall coeffi- 
cient RhNc/V, where N is the number of Cu atoms in a 
unit cell and V is the volume of the unit cell, can be used 
as a tool to identify the p valued Figure 1(b) compares 
the R H Ne/V of the Dy-doped Bi2212 samples (Ref. H) 
with those of LSCO (Ref. |43T ). from which the p values 
can be easily estimated with rather small uncertainties. 
Figure 1 (c) shows the obtained T c vs p relation for several 
pure and Dy-doped Bi2212 single crystals^ Note that 
the doping dependence of T c shows some deviations from 
the empirical formula T c /T™ ax = 1 - 82.6(p - 0.16) 2 , 45 
which is often used for determining the hole concentra- 
tion in Bi2212. In particular, the boundary between non- 
superconducting and superconducting phases turns out 
to be located at p w 0.07. 



B. Temperature dependence of thermal 
conductivity at "high" temperature 

Figure 2 shows the thermal conductivities of several 
pure and Dy-doped Bi2212 single crystals from 2 to 
150 K. The optimally-doped sample shows a well-known 
broad "QP peak" below T c as well as a "knee" feature 
around 10 K, whose origin has never been clarified. The 
size of QP peak (the enhancement from the minimum 
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near T c to the peak is 1.1 W/Km or 25 %) is among 
the largest in Bi2212 single crystals that have ever been 
studied , 17 i 22 ' 26 i 46 confirming the high quality of the crys- 
tals. In pure samples, either increasing or decreasing the 
hole concentration results in the weakening of the QP 
peak while the magnitude of total thermal conductiv- 
ity monotonously decreases with decreasing the hole con- 
centration. Furthermore, the peak suppression is much 
more significant in the underdoped samples than that in 
the overdoped samples, which is possibly related to the 
nanoscale inhomogeneity found by STM experiments in 
the superconducting state of this compound i 47 ' 48 ! 49 This 
behavior is very similar to that in the Zn-doped Bi2212, 22 
in which local suppression of the superconductivity lead- 
ing to a patchy superconducting state is playing a role^ 
Upon doping Dy, the QP peak is also strongly suppressed 
but a new peak at 10 K emerges. Interestingly, in the 
Dy81 (x = 0.17) sample, both the QP peak and the 10 
K peak show up, which means that the 10 K peak is not 
directly related to the QP heat transport. This low-T 
peak is the strongest in the Dy45 (x — 0.30) sample and 
survives in the non-superconducting sample (DyO, x — 
0.53), which clearly indicates that it originates from the 
phonon heat transport. 

In the undoped compound of other cuprates, like 
La2CuC>4, YBa2Cu306, Nd2CuC>4 and Pri. 3 Lao.7Cu04, 
the phonon heat transport behaves similarly to that 
of common insulators and shows a large phonon peak 
around 20 K, whose magnitude varies between ~ 20 
and ~ 80 W/Km.^i^^ However, it is known that 
Bi2212 system has much dirtier phonon heat transport 
than other cuprates^ 6 - probably due to the strong dis- 
order of the crystal lattice caused by the excess oxygen 
and the incommensurate superlattice along the b axis. 
The incommensurate modulation had been discussed to 
be related to the extra oxygen locating in BiO double 
layers and the lattice mismatch between the BiO block 
and the CuC>2 layer.— In an earlier work, the thermal 
conductivity of Bi2212 parent material was studied in 
Bi2Sr2YCu208 and it was found that the phonon heat 
transport is too weak to give rise to the phonon peak 
(only a hump-like feature observed at 20 - 30 K with 
the magnitude of k smaller than 2 W/Km). 46 To our 
knowledge, the present result shows the phonon peak 
of Bi2212 for the first time. It is not surprising that 
the phonon peak is located at much lower temperature 
than in other cuprates if the defect scattering is much 
stronger in Bi2212. The present result also indicates that 
a moderate Dy doping can enhance the phonon conduc- 
tivity rather strongly despite that the atomic substitution 
always introduces additional impurity/defect scattering 
on phonons. Apparently, the crystal structure is some- 
what stabilized by doping an appropriate amount of Dy. 
A possible scenario is that substituting Ca 2+ ions with 
Dy 3+ increases the oxygen content and results in the ex- 
pansion of the BiO bolck, which is helpful for relaxing 
the lattice mismatch between the BiO block and Cu02 
planed Note that the low-T peak observed in the Dy- 
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FIG. 3: (color online) Thermal conductivity of Dy-doped 
Bi2212 (x = 0.30) single crystal (annealed at 800 °C in air) 
and the same sample after oxygenating heat treatment (an- 
nealed at 500 °C in air). 

doped samples clarifies that the "knee" feature of k(T) 
in superconducting Bi2212 single crystals is due to the 
competition between the decrease of QP heat transport 
and the increase of phonon heat transport upon lower- 
ing temperature across 10 K. Furthermore, the phonon 
heat transport in pure Bi2212 is obviously very weak, and 
its magnitude at 10 K (where the phonon conductivity 
peaks) is smaller than or at most comparable to that in 
the non-superconducting Dy-Bi2212. 

Similar to the pure Bi2212 crystals, the heat trans- 
port of Dy-Bi2212 is also sensitive to the oxygen con- 
tent. After appropriate annealing to increase the oxy- 
gen content, for example, the x = 0.30 crystal (Dy45 in 
Fig. 2) presents a parallel-shift-like enhancement of the 
thermal conductivity at first glance, as shown in Fig. 3. 
However, we should note that the thermal conductivity 
enhancement above T c is mainly caused by the increase 
of the electron heat transport (which can be roughly es- 
timated from the normal-state resistivity data using the 
Wiedemann- Franz law), while the enhancement of the 
lOK-peak should be mainly due to the improvement of 
phonon transport since the QP (or electron) heat trans- 
port in this temperature region is expected to be strongly 
damped in highly Dy-doped samples. It is likely that the 
oxygenating process reduces the number of oxygen va- 
cancies in the crystal, which scatter phonons. The mag- 
nitude of the phonon peak at 10 K increases from 3.8 to 
5.2 W/Km, but both curves do not show the QP peak al- 
though the superconducting temperature increases from 
45 to 55 K. 



C. Temperature dependence of thermal 
conductivity at subkelvin temperature 

To further investigate the doping dependence of QP 
heat transport, the thermal conductivities of Bi2212 and 
Dy-doped Bi2212 single crystals are measured at very low 
temperatures down to 70 mK. Some of the data have been 
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reported in a previous publication (see Figs. 3 and 4 in 
Ref. 0). It has been shown that in the lowest-T regime 
where phonons are in the boundary scattering limit, the 
QP term and the phonon term of thermal conductivity 
can be separated by fitting data (below ~ 130 mK) to 



K 

T 



T 



bT 2 , 



(2) 



where the residual term kq/T is the QP contribution 
and bT 2 is the phonon contribution . 4 ' 9 ' 10 ! 14 Note that the 
residual thermal conductivity of a d-wave superconductor 
was discussed to present a "universal" behavior, that is, 
it is independent of the QP scattering rate if the impurity 
bandwidth satisfies ksT <C 7 <C Ao (Ao is the maximum 
gap)^ In this picture, the QP heat conductivity (n e /T) 
of cuprate superconductor decreases with temperature 
and finally arrives at a certain value, which depends only 
on the doping level, at low enough temperature. How- 
ever, our Bi2212 crystals clearly show much larger resid- 
ual thermal conductivity compared to the samples used 
in earlier works££ Judging from the methods of growing 
single crystals and the normal-state resistivity data^ it 
is obvious that the difference in kq/T is strongly related 
to the cleanliness of the crystals, which actually chal- 
lenges the validity of the "universal" behavior. Indeed, 
this issue was recently elucidated by a critical exami- 
nation of the "universal" behavior in a well-controlled 
experiment^ it was found that the disorder or inhomo- 
geneity of cuprates, which are especially strong in Bi2212, 
make the "universal" behavior to be only a rough approx- 
imation. Nonetheless, the residual thermal conductivity 
universally shows a monotonous doping dependence in 
all cuprates. 

Recently, a controversy arose for the data analysis of 
the very-low-T thermal conductivity of high-T c cuprates; 
namely, whether the commonly observed T-dependence 
of k/T data below 1 K weaker than T 2 is due to the spec- 
ular reflection of phonons off the sample surface or sim- 
ply because the phonon boundary scattering can only be 
achieved at temperatures lower than ~ 130 mK 4 ^ ' 11 ' 14 
In this regard, although Taillefer's group was the first 
to propose the analysis of Eq. @ and successfully de- 
scribed the very-low-T data of YBCO^ they have re- 
cently shifted to conjecture the existence of specular re- 
flections in a wide temperature range and analyzed their 
data of YBCO and LSCO by fitting to 



K Kn 



(3) 



where the second term with 7 < 2 is purported to be 
the phonon conductivity in the specular scattering limit J 
However, it was recently pointed out^ that one can es- 
timate the phonon mean-free path £ p in the temperature 
range of interest, which shows that that l v is too short 
for the specular-reflection scenario to be valid. Further- 
more, the analysis involving Eq. (3) analysis is actually 
based on another assumption that the QP thermal con- 
ductivity is independent of temperature up to at least ~ 
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FIG. 4: (color online) Low-T thermal conductivity of pure 
and several Dy-doped Bi2212 single crystals measured down 
to 70 mK, plotted as k/T vs T. The Dy contents x for Dy80, 
Dy42 and DyO samples are 0.20, 0.34 and 0.51, respectively. 
Strong downturn of the low-T data of OD83 sample is caused 
by the electron-phonon decoupling. 



1 K, which is inconsistent with an experimental finding 
from the same group^ that in YBCO the QP thermal 
conductivity at finite temperature increases quickly with 
temperature, following a T 3 dependence^ In this sense, 
the "good" fitting of Eq. Q to some experimental data 
are not physically meaningful. 

If one tries to analyze the data of Bi2212 using Eq. 
([3]) , it can be found that the low-T thermal conductivity 
actually follows 



- = — + AT, 
T T 



(4) 



below 0.6 K down to nearly 100 mK, as shown in Fig. 
4; here the power of T in the second term is clearly too 
small to be attributed to the phonon specular reflection of 
phonons^ Note that the same temperature dependence of 
K has recently been observed in Tl2Ba2CuOg+5 (T12201) 
system, whose doping levels are located in the overdoped 
regime. 8 Hawthorn et al. attributed the phonon conduc- 
tivity of AT to the dominant electron-phonon scatter- 
ing with high carrier concentration of T12201,— like the 
case in a metal, which actually means that the phonon 
boundary scattering limit can never be established in this 
system. Our present results, on the other hand, are ob- 
served mostly in the underdosed or optimally doped re- 
gions of Bi2212, indicating that the AT term may not 
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be the simple phonon contribution. Considering both 
the much larger k /T term of Bi2212 compared to other 
cuprates and the much weaker phonon transport seen in 
the "high"-T region (Fig. 2), it is reasonable to con- 
clude that the T 2 dependence of the low-T thermal con- 
ductivity (above ~ 100 mK) of Bi2212 is related to an 
amorphous-like phonon transport 5 ' 9 together with a In- 
dependent QP heat conductivity. This is essentially con- 
sistent with a previous finding in Bi2212 by Nakamae et 

To determine the residual term kq/T at zero temper- 
ature by employing some extrapolation of the data at 
finite temperature, one should always keep in mind that 
the formula used for the extrapolation must have a phys- 
ical reasoning, which assures that the same functional 
form holds throughout the range of the extrapolation (to 
K). Therefore, although Eq. (J) can fit the Bi2212 
data pretty well in a broad temperature range, this fit- 
ting cannot provide the precise determination of kq/T, 
because the T 2 -dependent QP heat conductivity does not 
hold at very low temperature (~ 100 mK or lower) when 
the temperature scale becomes smaller than the impurity 
bandwidthi££ For this reason, it is still the most reliable 
way to obtain kq/T by fitting the lowest-T data to Eq. 
(H]) , as we did in Ref. UA 

In passing, it should be pointed out that the very-low- 
T data of overdoped Bi2212 show a significant downturn 
below 0.2 K, which is known to be due to the electron- 
phonon decoupling at very low temperatures^ that is 
often significant in the overdoped cuprates^ This phe- 
nomenon further demonstrates that the T 2 behavior of 
low-T thermal conductivity in this system does not have 
the same origin as that in T12201, which was discussed 
to be the strong electron-phonon scattering^ 
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D. Magnetic-field dependence of thermal 
conductivity 

The magnetic-field dependence of both pure and Dy- 
doped Bi2212 crystals are studied in detail for H \\ c. 
Figures 5(a) and 5(b) show the results of OP95 sample 
at 5 - 100 K, which are essentially consistent with previ- 
ous studies i 17 ' 22 i 26 Note that the famous "plateau" fea- 
ture is absent in these k(H) isotherms since the thermal 
conductivity is always measured along the a axis in the 
present work. As shown in Fig. 6, the field-dependence 
of k can be well described by the simple relation 



K (H,T) 



k e (T) 



p(T)H 



Kph{T), 



(•5) 



where n e and K p h are the electronic thermal conductiv- 
ity and the phononic thermal conductivity, respectively; 
here the field dependence is attributed to the quasipar- 
ticle scattering off vortices^ Apparently, the amplitude 
of the field dependence of k is dependent on the contri- 
bution of QPs to the total thermal conductivity. One 
should notice that the magnetic field dependence of QP 



FIG. 5: (color online) Magnetic-field dependences of the low- 
T thermal conductivity of the pure Bi2212 single crystals, 
whose doping levels cover from underdoped to overdoped re- 
gions. The field direction is along the c axis. 



thermal conductivity in d-wave cuprates is determined by 
the competition between the decrease of heat conductiv- 
ity due to the vortex scattering and the increase of heat 
conductivity due to the increase of the nodal QP exci- 
tations, the so-called "Volovik" effect ™£L2Z&21£l The 
former effect is dominant at high temperature whereas 
the latter one is more important at low temperature. 
Thus, the weakliness of the field dependence of k at low 
temperatures could be due either to the decrease of QP 
portion of heat conductivity or to the "Volovik" effect, 
or both. 

As shown in Figs. 5(c-f), the field dependence of k is 
also rather strong in the overdoped sample, almost the 
same as that in the optimally-doped sample, whereas it 
is weakened by several times in the underdoped sample. 
Nevertheless, the k(H) data in these pure Bi2212 crystals 
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FIG. 6: (color online) Representative k(H) isotherms of pure 
Bi2212 single crystals with different doping level and the cor- 
responding theoretical fitting curves using Eq. ([5]). 



at all temperatures studied are well described by Eq. ([5|) 
(see Fig. 6 for several representative data), which means 
the QP scattering mechanism does not change upon tun- 
ing the oxygen content. Apparently, the weakening of 
k(H) in underdoped sample is due to the suppression of 
the QP contribution to the heat conductivity, consistent 
with the conclusion from the temperature-dependence 
data. 

Note that the above discussion is based on an assump- 
tion that the phonon thermal conductivity is indepen- 
dent of magnetic field, which however has recently been 
found to be not always valid for high-T c cuprates. In a 
parent compound of high-T c cuprates, PLCO, in which 
the phononic heat transport is very clean, diluted free 
spins associated with a small oxygen nonstoichiometry 
introduces strong magnetic-field dependence of phonon 
thermal conductivity through the phonon scattering off 
free spins Since the localized free spins are commonly 
found in superconducting cuprates, as is evidenced by the 
Schottky anomaly in the specific heat,^^^^^^. 
some fraction of the magnetic-field dependence of n is 
likely due to the spin-phonon scattering. In the case of 
pure Bi2212, however, the effect of spin scattering on the 
phonon heat transport would be unimportant because 
the phonons are strongly scattered by the lattice disor- 
der or defects. 

The magnetic-field dependences of the thermal con- 
ductivity of Dy-doped samples with x = 0.17, 0.30 and 
0.51 are shown in Fig. 7. Similar to the case of Dy- 
free crystals, the field dependence of k is gradually weak- 
ened with decreasing carrier concentration, which should 
be also due to the reduction of the electronic contribu- 
tion to the heat conductivity. An interesting point is 
that the k in Dy-Bi2212 samples show more complicated 
magnetic-field dependences. For x = 0.17, the magnetic- 
field dependence of n can still be well described by Eq. 
([5]). On the other hand, for x — 0.30, an "S"-shaped 
feature at 1 - 3 T shows up in the k(H) isotherms below 
30 K; the fitting to these data using Eq. ([5]), as shown 
in Fig. 8, suggests that there is an additional dip-like 
change of k superimposed on the background of usual 
field dependence described by Eq. (O and the position 



1.00 
0.98 
0.96 
0.94 
0.92 



1 1 1 I 1 I 


■ i 1 i ■ i 
— • — 5 K 




— • — 10 K 




— • — 1 2.5 K 


(a) \ 




- — i— 18.5 K 




24 K 




' — ^30K 




•— 40K 





1.00 


i i i i i i i i i i i i 


0.98 




90 K 
80 K_ 


0.96 






0.94 






0.92 


- (b) 


N. 50 K • 
40 K\ " 



1 2 3 4 5 6 
H(T) 



1 2 3 4 5 6 
H(T) 



1.00 



Dy0.30-Bi2212, T, = 45 K Dy0.51-Bi2212, NonSC 

1.00 



0.98 



•0.96 



0.94 




0.99 



0.98 



0.97 



0.96 




1 2 3 4 5 6 
H(T) 



FIG. 7: (color online) Magnetic-field dependences of the low- 
T thermal conductivity of the Dy-Bi2212 (x = 0.17, 0.30 and 
0.51) single crystals, which are superconducting at 81 K, 45 
K and non-superconducting, respectively. 
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([5]). The dashed line indicates the deviation of experimental 
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of the dip slightly moves to lower field with decreasing 
temperature. For the x = 0.51 sample, in which the over- 
all field dependence is much weaker, a step-like decrease 
of k clearly occurs around 2 T at low temperatures. 

To probe whether the peculiar k(H) behavior for 
higher Dy dopings is caused by the oxygen defects, the 
field dependences of k are also measured for the x = 
0.30 sample after re-annealed at 500 °C in air, which in- 
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FIG. 9: (color online) Magnetic-field dependences of the low- 
T thermal conductivity of x = 0.30 sample after re- annealed 
at 500 °C in air. 



creases the oxygen content and increases the T c from 45 
K to 55 K (the temperature dependence of k is shown in 
Fig. 3). It can be seen from Figs. 7(c) and 9 that the 
profiles of k(H) do not change so much with increasing 
oxygen content except that the dip-like feature becomes 
more pronounced at low temperatures. The position of 
the dip is still located at 1 - 3 T. This indicates that the 
oxygen defects are not directly related to the physical 
origin of this peculiar feature. 

Since both the quasiparticle transport and the oxygen 
defects are not likely to be producing the additional con- 
tribution to the field-dependence of k, it is natural to 
conclude that the Dy 3+ ions are playing the key role. 
Apparently, the x = 0.30 sample does not have any long- 
range order of Cu spins or Dy spins, so the Dy 3+ ions 
are most likely affecting the heat transport as free spins. 
A common magnetic scattering on phonons is caused by 
free spins or paramagnetic moments through spin-phonon 
interaction.^ In this case the energy splitting of the spin 
states that induces resonant scattering on phonons is in- 
creased by the Zeeman effect. The phonon scattering 
off these spins is most effective in suppressing the heat 
transport when the energy splitting is equal to ~ 3.8 
ksT, where the spectrum of phonon heat conductivity 
peaks j 16 i 59 therefore, the spin-phonon scattering usually 
generates a dip-like feature in k(H) at this energy and 
the dip position shifts to higher fields with increasing 
temperature^ As we have mentioned, the effect of the 
magnetic phonon scattering on the heat transport has 
been known for a long time and has recently been found 
to affect the field-dependence of k in cuprates.— While 
the magnetic scattering on phonons seems unimportant 
in pure Bi2212 single crystals, in which the phonon heat 
transport is strongly damped by defects and lattice dis- 
orders, it becomes much more pronounced upon doping 
Dy that not only brings free spins but also enhances the 
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FIG. 10: (color online) Magnetic-field dependences of the low- 
T thermal conductivity of the Dy-Bi2212 (x = 0.53) single 
crystal down to 0.36 K. The magnetic field is applied along 
the c axis. 



phonon heat transport. 

As shown in Fig. 7(d), a step-like k(H) isotherms is ob- 
served at low temperatures in the x = 0.51 sample. Is this 
also due to the phonon scattering by free spins? In this 
scenario, if the energy splitting is larger than 3.8 kgT, the 
magnetic field further increases the energy splitting and 
simply weakens the scattering and produces a step-like 
increase of thermal conductivity?^ Apparently, the step- 
like decrease of k(H) in Dy-Bi2212 is out of the expecta- 
tion from the common spin-phonon scattering mechanism 
and suggests some other physical mechanism. Since the 
x = 0.51 sample is non-superconducting with an insulat- 
ing resistivity behavior at low temperatures, the step-like 
k(H) must also be irrelevant to the electron transport. 



E. "Plateau" of low-temperature k(H) isotherms in 
Dy-Bi2212 

In Fig. 7(d), the k(H) isotherm of x — 0.51 sample 
at 5 K shows a step-like decrease at low field and n be- 
comes nearly independent of magnetic field above 2 T, 
a "plateau" -like behavior. We note that it is not clear 
whether this "plateau" behavior is the same as that ob- 
served in pure Bi2212 along the b axis^ because the 
latter is believed to have its origin in the QP transport, 
while in the present case the sample is an insulator. To 
investigate this intriguing behavior in Dy-doped Bi2212, 
the magnetic-field dependence of k are studied in more 
detail at very low temperature down to 0.36 K and in 
high magnetic field up to 16 T (along the c axis). As 
shown in Fig. 10, the "plateau" -like behavior is well 
reproduced in another sample of x = 0.53 at low tem- 
peratures. With lowering temperature, the step-like sup- 
pression of k becomes more significant and the strongest 
(~ 16%) suppression is observed between 0.97 and 0.5 K. 
At even lower temperature of 0.36 K, the field-induced 
suppression of k is getting weaker. The characteristic 
field, where the suppression completes or the "plateau" 
starts, is weakly dependent on the temperature, increases 
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FIG. 11: (color online) Temperature and field dependences 
of magnetization of a non-superconducting Dy-Bi2212 single 
crystal with x = 0.53. 



slightly from - 1 T at 0.36 K to - 2 T at 5 K. Somewhat 
intriguingly, this trend is the same as that found for the 
plateau in the superconducting Bi2212 samples^ 

Figure 11 shows the temperature and field dependences 
of magnetization for ani = 0.53 sample. It is clear that 
the magnetization data show rather conventional spin po- 
larization of paramagnetic spins, which must be spins of 
Dy 3+ ions. The characteristic field for spin polarization 
is about 2 T at 2 K and increases quickly with increasing 
temperature; this is rather different from the character- 
istic field of k(H) where the "plateau" appears. In this 
sense, the k(H) behavior does not appear to have direct 
relation with the Dy 3+ free spins and their polarization, 
which is different from the case of the x — 0.30 sample. 

It is known that in rare-earth- or Y-doped Bi2212, 62 
the antiferromagnetic (AF) order of Cu spins starts to 
be established when the superconductivity is completely 
suppressed and charge carriers are localized upon increas- 
ing the concentration of rare-earth ions. Therefore, one 
possible picture of the plateau behavior would be the 
magnon heat transport being affected by applied mag- 
netic field: In a long-range-ordered spin system, the low- 
energy magnon excitations can act as heat carriers, but 
when the magnetic field drives the magnon states to suf- 
ficiently high energies, these magnons are no longer ther- 
mally excited and drop from heat conduction^ 3 - This 
causes a plateau-like k(H ) behavior in magnetic systems, 
as has been observed in the yttrium iron garnet systemi^ 3 . 
However, one should keep in mind that the magnons are 
possibly involved only when the antiferromagnetic state 
has un-gapped (or very weakly gapped) acoustic magnon 
branches, which needs to be confirmed by neutron mea- 
surements. More importantly, the very abrupt nature of 
the onset of the plateau observed at very low temper- 
ature [Fig. 10(a)], together with the rather weak tem- 
perature dependence of the magnetic field from which 
the plateau starts, seems to be difficult to be understood 
within this picture. Also, it is even not clear whether the 
"plateau" is only observable in non-superconducting Dy- 



doped samples. In other word, although the "plateau" 
feature is not present in the low-T n(H) curves for the 
a-axis Bi2212 samples in strong fields up to 16 T^ the 
present data could not clarify whether it is also absent 
in high magnetic fields for superconducting (and not AF 
ordered) Dy-Bi2212 with lower concentration of x. De- 
tailed future study of this plateau phenomenon is clearly 
called for. 

In passing, we note that the "plateau" behavior grad- 
ually weakens with increasing temperature, which is pos- 
sibly due to the disappearance of the antiferromagnetic 
order. However, at even higher temperatures, the field 
dependence becomes again stronger above 18.5 K, as is 
shown in Fig. 7(d). The strongest ii-dependence is ob- 
served at 30 K, where the magnetic field of 6 T induces 
~ 3 % suppression of k without any signature of sat- 
uration. This magnetic-field dependence of k at high 
temperature is actually difficult to understand because if 
phonon transport is not affected by magnetic field, the 
only way to bring the field dependence of k is through the 
electron channel. In this regard, the magnetoresistance 
measurements have shown that the in-plane resistivity 
changes by less than 0.1 % in magnetic field (|| c) up 
to 14 T. 44 Such small change of resistivity cannot bring 
several percent change of thermal conductivity unless the 
Wiedemann-Franz law is strongly violated in this case. 



IV. SUMMARY 

Both the quasiparticle and phonon heat transport of 
Bi2Sr2CaCu20g+,5 are carefully studied for a broad dop- 
ing range by using pure Bi2212 single crystals with tuned 
oxygen contents and Bi2Sr2Cai_ a; Dy a ;Cu208+5 single 
crystals with different Dy contents. It is found that 
the QP transport decreases with lowering doping by ei- 
ther decreasing oxygen content or increasing Dy content, 
whereas the phonon conductivity is enhanced in mod- 
erately Dy-doped samples where a phonon peak is ob- 
served at 10 K, which means that Dy 3+ ions not only 
act as impurities or point defects but also stabilize the 
crystal lattice. The magnetic-field dependence of k in 
pure Bi2212 at temperature above 5 K is mainly due to 
the QP scattering off vortices. In contrast, the Dy-doped 
samples (with x = 0.30) present an additional "dip"-like 
term of k(H) at low field, which is discussed to be related 
to the phonon scattering by free spins of Dy 3+ ions. In 
non-superconducting Dy-Bi2212 samples with x ~ 0.50, 
a "plateau" -like k(H) isotherms with characteristic field 
at 1 - 2 T is observed. We discuss that this "plateau" is 
possibly related to the magnon transport. 
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